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ABSTRACT 



I study the possibility that a cooling flow (CF) exists at the main phase of 
super massive black hole (SMBH) growth during galaxy formation. To ensure 
that jets launched by the SMBH efficiently expel gas from the galaxy, as is 
required by recent results, the gas should be in the hot phase, rather than in cold 
clouds. The short radiative cooling time of the hot gas leads to the formation 
of a CF, but heating by the active galactic nucleus (AGN) prevents catastrophic 
cooling. Cold blobs that start as instabilities in the hot phase feed the SMBH 
from an extended region, form an accretion disk, and lead to the formation of 
jets. These jets can expel large quantities of gas out of the galaxy. This cycle, 
that is termed cold feedback mechanism in CFs in clusters of galaxies, might 
explain the correlation of SMBH to bulge masses. Stars are formed, but at a 
lower rate than what is expected when heating is not included. Such a CF is 
termed a moderate CF. 



1. INTRODUCTION 

The tight correlation between the supermassive black hole (SMBH) mass, Mbh, and the 
velocity dispersion, a, of the hot component of the host galaxy (e.g., Merritt & Ferrarese 
2001; Gebhardt et al. 2000; Giiltekin et al. 2009 ) is well established, and so is the correlation 
between the SMBH mass and the bulge mass, Mbui ge (e.g., Kormendy & Richstone 1995; Laor 
2001). It is possible that the M B h — M bu \ ge correlation is determined by a feedback process 
of the active galactic nucleus (AGN) of the host bulge, powered by the accreting SMBH. 
The feedback mechanism where AGN jets (outflow; wind) suppress gas from cooling to low 
temperatures and from forming stars was discussed for both cooling flows (CFs) in galaxies 
and clusters of galaxies (e.g., Binney & Tabor 1995; Nulsen & Fabian 2000; Reynolds et al. 
2002; Omma & Binney 2004; Soker & Pizzolato 2005), and in galaxy formation (e.g., Silk 
& Rees 1998; Fabian 1999; King 2003; Croton et al. 2006; Bower et al. 2008; Shabala & 



1 Dept. of Physics, Technion, Haifa 32000, Israel; soker@physics.technion.ac.il 



-2- 



Alexander 2009). By galaxy formation I will refer also to bulge formation in a spiral galaxy, 
and vice versa. 

There are other attempts to explain this correlation. For example, in their model of 
bulge formation by coalescence of giant (M ~ 1O 8 M ) clumps, Elmegreen et al. (2008b) 
considered feedback by supernovae (SNe), but not by AGN activity. In their model, the 
SMBH at the center of the galaxy is formed from the merger of the intermediate mass BHs 
that reside in each clump (Elmegreen et al. 2008a). The ratio of SMBH to bulge mass 
M B H/^fbui g e ~ originates in the massive clumps. In the present paper I assume that 
the SMBH forms by accreting mass from the ISM, and that the Mbh — ^tmige correlation is 
determined by an AGN feedback mechanism. It is determined in the sense that the AGN 
activity limits the ISM mass that is eventually converted to stars by expelling it out of the 
galaxy. The ejection of large quantities of gas out of the galaxy during its formation was 
found to be a necessary process in recent studies (e.g., Bower et al. 2008). 

In a previous paper I tried to account for the SMBH— bulge masses correlation with 
a feedback mechanism based on jets launched by the SMBH. This feedback is based on 
narrow jets that are launched by the central SMBH, and expel large amounts of mass to 
large distances. The condition for an efficient expelling process is that the jets do not 
penetrate through the inflowing gas, such that they can deposit their energy in the inner 
region where the bulge is formed. A relation between the mass accreted by the SMBH and 
the mass that is not expelled, and is assumed to form the bulge, was derived (Soker 2009; 
this derivation is repeated in section 4.3). It was noted that the same mechanism could 
operate in suppressing star formation in CF clusters, making a tight connection between the 
feedback in galaxy formation and CF clusters. In the present paper I extend the comparison 
between the feedback mechanism operating in CF clusters and during galaxy formation. My 
fundamental assumption is that the correlation between the SMBH mass and the host galaxy 
bulge mass is determined by an AGN feedback process that expel large quantities of gas out 
from the galaxy. 

In section [2] and [3] I argue that for the feedback mechanism to be efficient, as required 
from my fundamental assumption, most of the gas should pass through the hot phase, unlike 
the case in the model studied by Binney (2004, 2005). In section H] I discuss a model where 
the feedback between the SMBH and star formation in the bulge occurs through a process 
similar to that operating in CFs in clusters of galaxies. A comparison of CF in clusters to 
the process at galaxy formation was done in the past (e.g., Croton et al. 2006). In the 
model discussed by Binney (2004; 2005), for example, heating by the AGN activity and the 
ejection of large quantities of gas out from the galaxy are crucial processes, as they are in 
the present model. However, Binney (2004) proposed a model where the gas that forms star 
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has never been heated to the virial temperature, while that gas that was heated to the virial 
temperature has negligible contribution to star formation. In the model proposed in section 
14. 2^ on the other hand, large quantities of gas in the hot phase cool to form stars. The 
advantage in the presence of large amount of mass in the hot phase is that the hot phase 
is more susceptible to AGN heating (Binney 2004; Dekel & Birnboim 2006; Cattaneo et al. 
2006; Hopkins & Elvis 2009). I summarize the proposed model in section [5j 



2. PROBLEMS WITH FEEDBACK DRIVEN BY BONDI ACCRETION 

I define the Bondi accretion radius by 

2GM BU 6GM BU [im H 
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where Mbh is the SMBH mass , C s = (hkT ' j'in.mii) 1 ! 2 is the ISM sound sped, and all 
other symbols have their usual meaning. The Bondi accretion rate for an adiabatic index of 
7 = 5/3 is (Bondi 1952) 

M B ondi = 0.25nR 2 B pCs, (2) 

where p is the density at a large distance from the accreting body. The characteristic inflow 
time is 

— Rb (q\ 

Tflow = -pr- [o) 

(The inflow velocity at Rb is <C C s , and the practical inflow time scale is longer than r flow .) 

A condition for the existence of a Bondi type accretion flow is that the radiative cooling 
time is longer than the inflow time r coo i > Tfl ow . Otherwise the ISM rapidly cools, the value 
of C s rapidly decreases, and the assumptions that lead to the Bondi accretion flow in a hot 
ISM break down. Namely, we are in the regime of accreting cold gas, the problems of which 
are discussed in section [3j This condition reads 

3 nkT Rb 

2An e n p > ~C~ S ' ^ ' 

where n, n e , and n p are the total, electron, and proton number density, respectively. Condi- 
tion (j3J) can be cast in the form 

52 ^L RbC 2 > 0.257^0,. (5) 

For the cooling function A I assume a zero metallicity composition and take the temper- 
ature range T > 10 7 K. This gives A (T) ~ 2 x 10" 23 (T/10 7 K) 1 / 2 erg cm 3 s" 1 (Sutherland & 



-4- 



Dopita 1993). For lower temperatures and higher metallicity the conclusions will be stronger 
even, as cooling time is shorter. Substituting equation (jSJ) into equation (J2J) with the above 
approximation for A, gives 

Wk) J^ Mq yr_1 - (6) 

The actual limit is lower even for three reasons: (1) Some fraction of metals will be presence, 
increasing the cooling rate. (2) The inflow time should be shorter than the cooling time not 
only from the Bondi radius, but also from a distance of several Bondi radii. (3) In galaxies 
that do not sit at the center of groups and clusters, the virial temperature of the gas is 
< 10 7 K. For a temperature of T = 2 x 10 6 K and zero metallicity, the cooling function is 
lower by a factor of 1.6, and the cooling time shorter by a factor of 5/1.6 = 3. The numerical 
coefficient in equation flBJ) is ~ 5 x 1O~ 1O M yr -1 . 

We can compare the limit from the Bondi accretion as given in equation (JS)) with the 
Eddington limit 

M Edd = 2.5 x 1(T 8 (^) ^M. yr" 1 , (7) 

where e is the efficiency of converting mass to radiation in the accretion process. We find 
that Meondi < O.lMEdd- It seems there is no time for the BH to grow if it is limited by the 
Bondi accretion, as the e-folding time is 

Mbr « / T \ ~ 1/2 

T > = 5^ > 6 X 10 {Wk) - < 8 > 

To grow by a factor of 1000, the minimum time required is ~ 4 Gyr, with a more likely value 
of > 5 Gyr. Clearly, this process cannot explain the rapid growth of SMBH during galaxy 
formation. 

There are two possible solution to the feeding problem. 

(1) Most of the ISM mass is in a cold T < 10 4 K <C T viria i phase, for which the Bondi 
accretion radius is much higher. The problem with cold clouds is that they are very dense. 
It is very hard to remove such clouds from the galaxy. More over, cold clouds are more likely 
to form stars. This solution makes it very hard to explain the termination of star formation. 
This is the subject of the next section. In any case, models that include Bondi accretion 
consider the hot T > 10 6 K phase. 

(2) Most of the mass is in the hot phase with T ~ T virial . However, the SMBH is fead by cold 
clumps that are embedded in the hot ISM. The clumps as a group contain a small fraction 
of the total ISM mass. This solution to maintain feedback heating in cluster CFs is termed 
the cold feedback mechanism (Pizzolato & Soker 2005, 2010; Soker 2006, 2008b). In section 
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|4] I compare this mechanism as it operates in galaxy formation at high redshifts to its role 
in low redshift cluster CFs. 



3. PROBLEMS IN MAINTAINING FEEDBACK WITH MASSIVE COLD 

ISM 

It is possible that most of the ISM mass resides in cold (~ 10 4 K) clouds, rather 
than being close to the virial temperature and in a hydrostatic equilibrium. The problem 
in this case is that the AGN feedback efficiency is too low, because it is extremely difficult 
to expel dense clouds from the inner regions of the galaxy. In a recent attempt Hopkins & 
Elvis (2009) note that AGN feedback works more efficiently on the hot phase, and consider a 
chain of processes to overcome this problem. In their model the clouds contain ~ 90% of the 
mass, while ~ 10% of the ISM mass is in the hot phase at hydrostatic equilibrium. A shock 
propagating through the hot phase causes the clouds to expand, such that radiation pressure 
from the AGN is more efficient in expelling the bloated clouds. Even in this more efficient 
scenario, the mass expelled in their model is not sufficient for the feedback mechanism studied 
here. Hopkins & Elvis (2009) take the bulge mass to be 1O 11 M , while the ISM mass is only 
1O 1O M . After a time of t = 10 8 yr from the start of their calculation, the mass expelled 
in their model is AM e = 5 x 1O 9 M . During the same period, for an efficiency of 10%, the 
SMBH mass has grown by AM SMB h ^ 1O 8 M . This ratio of AM e /AM SMBH ~ 50, is lower 
by more than an order of magnitude from the one required in feedback models that account 
for the bulge to SMBH masses ratio, if the feedback is to determine the final mass of the 
bulge. Had they include radiative cooling, the efficiency of expelling the ISM would be lower 
even. 

The conclusion is that for an AGN feedback to work, most of the gas in the inner regions, 
r < 10 kpc ~ 0.1R V , must reside in the hot phase; _R V is the virial radius. 

Although the AGN feedback requires the inner regions to be mostly in the hot phase, 
the gas feeding the galaxy at larger regions can be cold. In the cosmological cold streams 
model (Dekel et al. 2009a, b and references therein) three cold, low entropy, streams of gas 
penetrate the virial radius without being shocked, and reach the central region. Shocks near 
the center are not resolved in the numerical simulations presented by Dekel et al. (2009a), 
as their resolution is 1.4 kpc. Eventually, the streams must encounter a shock wave near the 
center. The typical influx rate considered by Dekel et al. (2009a) is m ~ 3OM yr _1 rad 2 . 
At a distance r from the center the post-shock total (protons, electrons and nuclei) number 
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density is 

,3OM yr-^adV V300 km s_1 ^ V 5 k P c , 
where v is the preshock inflow velocity, scaled with the inflow velocity at r ~ 5 — 10 kpc 
from the center according to Dekel et al. (2009a). The postshock temperature is ~ 10 6 K. 

For a temperature of ~ 10 6 K, and zero metallicity, the cooling function is A — 
10~ 23 erg cm 3 s" 1 (Sutherland & Dopita 1993). The cooling time of the post-shock gas 
is 

r ~ 3 x 10 6 ( ™ ^ " ( V Y /2 (^-] " ( A ° ^ vr 

V3OM yrtadV UOOkms- 1 / \5 kpc J VlO" 23 erg cm 3 s"V ' 

(10) 

where the dependence of the post shock temperature on velocity has been used. The half 
opening angle of a stream is ~ 10 — 15°; I scale with 12.5°. The pressure of the post-shock 
gas is larger than that of its surroundings, and it expands to the sides in a typical time of 

r sin 12.5° R / v \ - 1 / r \ , , 

Tf- = 3.5xl0 6 ( r) — yr. 11 

' v V300 km s' 1 / \5 kpc J J K J 

At r ~ 5 kpc the radiative cooling time is about equal to the expansion time of the 
postshock region t cs ~ Tf. Heating by the AGN will make the formation of a shock wave more 
likely. Cantalupo (2010) showed that ionization by star forming regions can remove cooling 
agents (ions) from the gas and by that efficiently reduces radiative cooling and transforms a 
cold mode accretion into a hot one. The effect is larger even when radiation from the AGN 
is considered. The effect of removing cooling ions comes in addition to the heating discussed 
next. Let the gas mass in the inner region be about equal to the stellar mass ~ 1O 11 M . 
With a total number density of n ~ 1 cm" 3 the mass resides within a radius of ~ 10 kpc. 
To maintain this mass that has a radiative cooling time of ~ 3 x 10 6 yr at a temperature of 
T ~ 10 6 K, requires a heating power of En ~ 4 x 10 44 erg s -1 . With an efficiency of 10% in 
converting mass to energy, the accretion rate onto the SMBH is ~ 0.1 M yr" 1 . In 10 9 yr the 
SMBH gains a mass of ~ 1O 8 M . In reality, the amount of ISM mass will be lower, as part 
of it is continuously forming stars, while some of the incoming gas is expelled back to large 
distances. This crude calculation and the reduction in radiative cooling rate as discussed 
by Cantalupo (2010) show that heating by the AGN can maintain a shock wave at a radius 
of r > 5 kpc. 

There is no problem in the scenario proposed here that the streams penetrate the virial 
radius, at R v ~ 50 — 100 kpc, as in the simulations presented by Dekel et al. (2009a). 
However, the model here does require that the streams eventually are shocked at r > O.lRy, 
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and a hot pseudo static atmosphere is formed. The presence of hot static atmosphere around 
an already grown SMBH was assumed before, e.g., Short & Thomas (2009). 

4. A COOLING FLOW PHASE 

4.1. The cold feedback mechanism in low redshift clusters 

In the cold-feedback model for clusters of galaxies (Pizzolato & Soker 2005, 2010; Soker 
2006; 2008a) mass accreted by the central SMBH originates in non-linear over-dense blobs 
of gas residing in an extended region of r ~ 5 — 50 kpc; these blobs are originally hot, but 
then cool faster than their environment and sink toward the center (see also Revaz et al. 
2008). The mass accretion rate by the central black hole is determined by the cooling time 
of the ICM, the entropy profile, and the presence of inhomogeneities (Soker 2006). Most 
important, the ICM entropy profile must be shallow for the blobs to reach the center as cold 
blobs. Wilman et al. (2008) suggest that the behavior and properties of the cold clumps 
they observe in the cluster A1664 support the cold feedback mechanism. 

The cold feedback mechanism in clusters has the following consequences. 

(1) Cooling flows do exist, but at moderate mass cooling rates: the moderate CF model 
(Soker et al. 2001). Indeed, in many CF clusters the heating cannot completely offset cool- 
ing (e.g., Clarke et al. 2004; Hicks & Mushotzky 2005; Bregman et al. 2006; Salome et al. 
2008; Hudson et al. 2010) and some gas cools to low temperatures and flows inward (e.g., 
Peterson & Fabian 2006) and forms stars (Wise et al. 2004; McNamara et al. 2004). Star 
formation is prominent when the radiative cooling time of the hot gas is short (Rafferty et 
al. 2008). These observations suggest that indeed cooling of gas from the hot phase to low 
temperatures does take place, including star formation. It is important to note that star 
formation and the feeding of the SMBH occur simultaneously. 

(2) The cold feedback mechanism explains why real clusters depart from an 'ideal' feedback 
loop that is 100% efficient in suppressing cooling and star formation. Simply, the feedback 
requires that non-negligible quantities of mass cool to low temperatures. Part of the mass 
falls to small radii. Part of this mass forms star, another part is ejected back, and a small 
fraction is accreted by the central SMBH. 

(3) Part (likely most) of the inflowing cold gas is ejected back from the very inner region. 
This is done by the original jets blown by the SMBH (Soker 2008b). The ejection of this gas 
is done in a slow massive wide (SMW) bipolar outflow, which are actually two jets. The 
basic mechanism is that the jet does not puncture a hole in the ICM, but rather deposits 
its energy in the inner region. Wide enough jets deposit their energy in the inner regions 
rather than puncturing a hole and expanding to large distances. Rapidly precessing jets or a 
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relative motion of the medium also prevent such a puncturing. In the case of the formation 
of low density bubbles in the ICM, it has been shown that in addition to SMW jets, precess- 
ing jets (Sternberg & Soker 2008a; Falceta-Goncalves et al. 2010) and a relative motion of 
the medium (Morsony et al. 2010) also lead to the inflation of bubbles close to the center. 
The most striking example of a SMW outflow is presented in the seminal work of Moe et al. 
(2009; also Dunn et al. 2009). By conducting a thorough analysis, Moe et al. (2009) find 
the outflow from the quasar SDSS J0838+2955 to have a velocity of ~ 5000 km s~ x , and a 
mass outflow rate of ~ 6OOM yr _1 , assuming a covering fraction of 5 ~ 0.2. The cooling 
and ejection back to the ISM of large quantities of mass, make the feedback process not only 
of energy (heating), but also of mass (Pizzolato & Soker 2005). 

(4) Such SMW jets can inflate the 'fat' bubbles that are observed in many CFs, in clusters, 
groups of galaxies, and in elliptical galaxies (Sternberg et al. 2007; Sternberg & Soker 2008a, 
b). The same holds for precessing jets (Sternberg & Soker 2008a; Falceta-Goncalves et al. 
2010) or a relative motion of the medium (Morsony et al. 2010). 

(5) The same mechanism that form SMW jets in CF clusters and galaxies, can expel large 
quantities of mass during galaxy formation, and might explain the SMBH-bulge mass corre- 
lation (Soker 2009; section 4.3 here). 

Following the problems with maintaining AGN feedback at galaxy formation by the 
Bondi accretion (sec. |2]), and the discussion in section [31 I turn to elaborate on the possibility 
that the feedback at galaxy formation is similar to that in CF clusters. 

4.2. A moderate cooling flow (CF) phase at galaxy formation 

In their study of galaxy formation with AGN feedback, Croton et al. (2006) consider 
the formation of an extended (up to the virial radius of the dark halo) hot atmosphere, and 
that gas from the central region of this atmosphere might be accreted onto a central object 
through a CF. In this section I take upon this, and consider the moderate CF model with 
a cold feedback. Namely, the gas that feeds the SMBH originates in an extended region as 
cooling blobs, and not as a Bondi type accretion. Heating of the inner regions by the AGN is 
a significant process, and the AGN activity level is regulated by the accretion of cold blobs 
from an extended region. 

I differ from Croton et al. (2006) in a crucial manner. In their model the main growth 
of the SMBH occurs in the 'quasar mode' at high redshifts (z > 2). During the CF phase of 
Croton et al. (2006) the accretion rate is low, and the mass of the SMBH does not increase 
by much. As shown in section 4.3, in the present model, the correlation between the SMBH 
mass and the bulge mass is set by a feedback process. For the feedback process to be efficient, 
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I require the presence of a CF. Therefore, in the present model the CF appears much earlier 
than in the model discussed by Croton et al. (2006). The redshift at which the moderate 
CF operates to form stars and feed the SMBH is determined by the specific model of galaxy 
growth that one uses. I do not argue for a different global feeding of gas to the galaxy, and 
therefor the SMBH growth period, as well as that of the bulge, is as in the specific model, 
e.g., at z > 2 in the model of Croton et al. (2006), and in the redshift range of z ~ 1 — 5 in 
the model of Bower et al. (2006, 2008). I limit myself to claim that a CF phase takes place 
during the main growth phase. The proposed process can be incorporated in large scale 
simulations as sub-grid physics. Namely, when mass is flowing to feed the central region, for 
both star formation and SMBH growth, one can use the proposed mechanism to operate the 
feedback process (as the inner region is not resolved by the large scale simulations). 

I start when the SMBH mass is already ~ 0.1—10% of its final mass, and the bulge/galaxy 
is already in the formation process. However, most of the stars in the bulge (or elliptical 
galaxy) have yet to be formed, and so does the SMBH mass. Cold streams that feed the 
galaxy (Dekel et al. 2009a, b) are assumed to be shocked at r > 10 kpc, such that most of 
the mass resides in the hot phase. The cooling time of the hot phase is short, and thermal 
instabilities lead to the formation of cool blobs that fall inward. These cool blobs are not 
much cooler than their surrounding for most of their journey inward (Pizzolato & Soker 
2010). Namely, they are not as cold as the gas discussed in section [3) The cool blobs are 
not much denser than their hot surroundings, and they can be expel relatively efficiently by 
the jets launched by the AGN when AGN activity is high (see section 4.3). A fraction of 
these blobs feed the SMBH. Heating by the SMBH activity, mainly by jets, facilitated the 
formation of this structure. The jets launched by the accreting SMBH not only heat the 
gas, but as is the case in the cold feedback mechanism (Pizzolato & Soker 2005), the jets 
accelerate large quantities of gas outward. A structure similar to that in cluster CF at low 
redshifts has been formed. Table 1 compares the properties of the two types of the moderate 
CF models. 

There are, however, two prominent qualitative differences between the proposed CF 
model at galaxy formation and that in CF in clusters of galaxies. 

(1) In the case of galaxy formation the cooling time in the inner ~ 1 kpc is shorter than the 
inflow time (section^. The gas cools very rapidly, pressure support is lost, and an inward 
supersonic flow is formed (Soker & Sarazin 1988). In clusters, on the other hand, the cooling 
time is longer than the dynamical time at all radii, and no such flow is formed; the gas 
feeding the SMBH is cold, but it originates in the hot phase. Some discussion of this feeding 
mode is given by Croton et al. (2006). In addition to this supersonic flow, cold blobs formed 
from the hot phase at larger radii feed the SMBH as well, as in the cold feedback model for 
heating CF in clusters (Pizzolato & Soker 2005, 2010; Soker 2006; 2008a). 
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Table 1: Comparing Moderate Cooling Flow Models'- 1 -* 



Property^ 2,3 ) 


Low z Clusters 


Galaxy formation 


Central e~ density (n ce ) 


0.1 cm~ 3 


10 cm -3 


Central temperature (T c ) 


3 x 10 7 K 


10 6 K 


System age (T age ) 


10 10 yr 


10 9 yr 


Central cooling time^ 4 ) (r c ) 


5 x 10 8 yr 


10 5 yr 


Cooling radius^ 5 ** (r c ) 


100 kpc 


30 kpc 


Dynamical time at 






0.01r c (r dl ) 


10 6 yr 


5 x 10 5 yr 


T c /r rf i 


500 


0.1 


Raw cooling rate^ 6 ) (Mj) 


10 - 1O 3 M yr" 1 


10 - 100M o yr" 1 


Star formation rate (M*) 


0-O.lMa, 


~ 0.1 - 0.5M,, 


Source of gas feeding 


Cold blobs from an 


Cold blobs from an 


the SMBH 


extended region. 


extended region 






+ inner supersonic inflow. 


Fate of most gas expelled 


Inflating large 


Expelled from the galaxy. 


from the inner region 


low-density bubbles 






(a) Depart from 100% 


(a) Most of the ISM is 




efficiency in suppressing 


susceptible to SMBH jets. 


Results and implications 


star formation. 


(b) Expelling huge amounts 


of the cold feedback 


(b) Shallow entropy profile. 


of mass from the galaxy. 


mechanism 


(c) Massive outflows that 


(c) Can account for the 




can inflate 'fat bubbles'. 


Mbh — Mb u ige correlation 




(d) Can operate without 


(Soker 2009; §4.3). 




the need for failed 


(d) Can operate without 




Bondi accretion (^2]). 


the need for failed 






Bondi accretion. 



Table 2: (1) In the Moderate CF model heating is important, but cooling of gas to low temperatures does 
occur, although at a much lower rate than that expected if no heating exists. 

(2) The values are crude, and most are given to an order of magnitude. 

(3) Some of the time scales are calculated by using equations from section [2] 

(4) The cooling time in galaxy formation is lower even, as a zero metallicity was assumed here, but some 
metals will be present at an age of ~ 10 9 yr. 

(5) The cooling radius r c is the radius at which the radiative cooling time (no heating included) equals the 
age of the system. 

(6) Raw cooling rate is the mass cooling rate if no heating was presence. 
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(2) In clusters, the region where cooling takes place is ~ 10 — 50 kpc, or about a fraction 
~ 0.1 — 0.3 of the cooling radius. Most of the ICM mass resides outside this radius. Any 
AGN activity can move mass from inner regions to regions further out, but the huge amount 
of mass in the outer regions prevents the mass to flow to very large distances from the center. 
The situation is different in the proposed CF model at galaxy formation, where most of the 
mass reside in regions having a cooling time shorter than the age of the system. The CF is 
a major process. Still, the heating by the central AGN ensures the presence of a hot phase, 
and prolongs the time it takes the gas to cool. That most of the mass is residing in the hot 
phase ensures an efficient coupling of jets launched by the AGN to the ISM. This coupling 
will be particularly efficient when the mass inflow rate is > 10 3 times the accretion to the 
SMBH, a ratio that might account for the SMBH— bulge masses correlation (Soker 2009; 
section 4.3 below). 

It is interesting to compare the model proposed here with that discussed by Binney 
(2004). In both models heating, aided by reduction in radiative cooling rate (Cantalupo 
2010), by the AGN activity is crucial, and in both models the ejection of large quantities 
of gas out from the galaxy takes place. The main difference is that Binney (2004) attributes 
all star formation to gas that was never heated to about the virial temperature. Instead, I 
suggest that a large fraction of the gas is shocked to the virial temperature, and is further 
heated by the AGN to prolong its hot phase. Most of the gas stays in the hot phase for a 
time longer than the dynamical time, and forms a (pseudo) static medium. Large quantities 
of this gas later cool, and form stars, as observed in cluster CF (but at lower efficiency). As 
the hot gas is much more susceptible to AGN activity, it allows for a feedback process from 
the AGN to work, and determine the masses ratio of the SMBH and bulge (Soker 2009). 



4.3. Correlation of SMBH-bulge masses 

In this section I briefly summarize the derivation of the correlation between the SMBH 
mass and the host galaxy bulge mass (Soker 2009). I will not repeat all steps and will not 
explain all assumptions, as they are in that paper. I will, however, make some modifications 
to account for the present proposal of a CF phase at galaxy formation, and to incorporate 
the new results of Soker & Meiron (2010). 

The basic assumptions are as follows. 

(1) The feedback mechanism is driven by jets. 

(2) The properties of jets launched by SMBH have some universal properties. As is shown 
in a new paper (Soker & Meiron 2010), the basic property is the the ratio of the total 
momentum discharge in the jets to the quantity M acc c: e p = Mjt>j/(M acc c), where M acc is 
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the accretion rate onto the SMBH, Mf is the mass flow rate into the two jets, and Vf is the 
jets' speed. Soker & Meiron (2010) find from statistical analysis of tens of galaxies that in 
the feedback model for the correlation e p = 0.038 ± 0.06. 

(3) The mass flowing in at early stages, i.e., the mass available for star formation, is very 
large. Namely, the mass that is converted to stars is limited by the feedback mechanism and 
not by the mass available in the SMBH surroundings. This is supported by studies of galaxy 
formation (e.g., Bower et al. 2008). 

(4) There is a relative transverse (not to be confused with the radial inward velocity) motion 
between the SMBH and the inflowing mass of t> re i — o. 

(5) The cooling surrounding mass M s that resides at a typical distance r s and having a 
density p s (see below), is flowing inward at a velocity of ~ a. Thus, M s ~ 4nr^ap s , and 
it is resupplied on a time scale of ~ r s /a. This mass will mainly form stars if it is not 
expelled by the jets. This assumption is in the heart of the proposed CF phase at galaxy 
formation. Namely, that the inflowing gas, be it a cooling flow that reaches a fast speed at 
r ~ 1 kpc (where the cooling time is short) or composed of cool clumps, is vulnerable to 
the jets launched by the central SMBH. The inflowing gas or clumps, are cooler than the 
virial temperature but not by much, hence their density is not much lower than that of the 
surrounding hot gas. They are not the very dense clumps at low temperatures that were 
discussed in section [31 and they can be expel by the jets. 

We note that M s S> M acc , as only a small fraction of the inflowing gas at scales of 
~ 0.1 - 10 kpc is accreted by the SMBH. 

If the jets penetrate through the surrounding gas they will be collimated by that gas, 
and two narrow collimated fast jets will be formed, similar to the flow structure in the 
simulations of Sutherland & Bicknell (2007). By fast it is understood that the jet's velocity 
is not much below its original velocity. If, on the other hand, the jets cannot penetrate the 
surrounding gas they will accelerate the surrounding gas and form SMW (slow-massive- wide) 
outflow (Soker 2008). 

I now derive (Soker 2009) the conditions for the jets not to penetrate the surrounding 
gas, but rather form a SMW outflow. Let the jets from the inner disk zone have a mass 
outflow rate in both directions of Mf, a velocity Vf, and let the two jets cover a solid angle 
of 47r<5 (on both sides of the disk together). The density of the outflow at radius r is 

Mf , N 

3 AnSr 2 Vf 

Let the jets encounter the surrounding gas residing within a distance r s and having a typical 
density p s ; this is the inflowing cooling gas, that if it is not expelled will form stars. The 
head of each jet proceeds at a speed Vh given by the balance of pressures on its two sides. 
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Assuming supersonic motion this equality reads p s v\ = p/(f/ — Vh) 2 , which can be solved for 

v h 

\ 1/2 / . \ 1/2 

Aix5r 2 s Vfp, 
Mf 

1225 (^) ( ±X /2 ( V -l) 1/2 ( a —Y'\ (13) 

1 1 1 0.1 y \cJ V200 kms- 1 / V ; 

where in the second equality the mass inflow rate M s ~ 4irp s ar 2 (by assumption 5), has 
been substituted. The time required for the jets to cross the surrounding gas and break out 
of it is given by 

—. — - = 4 x 10 — yr, (14) 

M f a ) \l kpcj K 1 

where in the last equality the same values as in equation (Tl3|) have been used. 

If there are no changes in the relative geometry of the SMBH and inflowing mass, the 
jets will rapidly penetrate the surrounding gas and expand to large distances. In this case 
the jets will not deposit their energy in the inflowing gas. For an efficient deposition of 
energy to the inflowing gas, we require that there will be a relative transverse (azimuthal) 
motion between the SMBH and the inflowing gas, such that the jets continuously encounter 
fresh mass. The relevant time is the time that the transverse motion of the jet crosses it 
width r s = Dj/v TC \ ~ Dj/a, as by our assumption 4 the relative velocity is v TC \ ~ a. The 
width of the jet at a distance r s from its source is Dj = 2r s sin a, where a is the half opening 
angle of the jet. For a narrow jet sin a ~ a ~ (25) 1 / 2 , and 

v rel \1 kpcj V200 km s~ V \0.1 J 

The demand for efficient energy deposition, t s < t p , reads then 



wr 8 i- (16) 

This result can be understood as follows. The ratio Vf<j/v 2 el comes from the ratio of the 
ram pressure of the narrow jet to that of the ambient gas which disturbs the jet, and from 
the relative transverse motion of the jet and the abient gas. The number 8 comes from the 
geometry of a narrow jet with a relative transverse velocity to that of the ambient gas. Using 
the definition e p = M/t>//(M acc c) from assumption 2, we derive 



M s ^ „ crc inn / f„ \ / a \ _1 / a \ 



M acc ~ ~~ P v 2 rel V0.04/ V200 km s~V \v iei J 



2 



>8e p ^ = A80^ ——-- T — . (17) 



-14- 



Again, it is expected that in its formation phase the galaxy will not be fully relaxed, and 
that the relative transverse velocity of the AGN and the inflowing gas will be of the order 
of the stellar dispersion velocity, i.e., t> re i ~ a 

The accretion rate M acc is the accretion rate onto the SMBH, and the inflow rate of the 
surrounding gas is assumed to form stars in the bulge (if it is not expelled by the jets). If 
the inflow rate is above the value given by equation ( ITT)) , the deposition of energy by the 
jets is efficient enough to expel the mass back to large distances and heat it (Soker 2008b). 
The interaction of the (narrow or wide) jets blown by the SMBH with the inflowing gas will 
form a wide outflow (SMW jets), that will expel more of the hot gas that is vulnerable to 
the jets. Namely, the jets blown by the SMBH will not allow the bulge to form stars at a 
rate larger than the value of M s given by equation ( ITT)) . Following Soker (2009) then, the 
SMBH to bulge mass ratio is about equal to M acc /M s . Equation ( ITT)) yields 

M BH ^ 0.002M bulEC ( -^-Y 1 ( ^ -) (—\ . (18) 

g V0.04/ V200 km s- 1 / V^rei/ 

The last equation closes the feedback cycle, in showing that a correlation can be driven by 
jets blown by the SMBH into the hot ISM. A key issue is that the medium is in the hot 
phase such that its density is not too high, and therefore it is vulnerable to the action of the 
jets. This hot phase feeds the SMBH via the process of a moderate cooling flow. 



5. SUMMARY 

Results from recent years show that the process of galaxy formation requires non- 
gravitational energy source not only to heat the gas, but also to expel large quantities 
of gas out from the galaxy (Bower et al. 2008). To efficiently eject the ISM from the galaxy 
by AGN activity the gas must be in the hot phase, namely, its temperature must be about 
the virial temperature (e.g., Hopkins & Elvis 2009). The conclusion from these studies is 
that most of the ISM during galaxy formation must evolves through the hot phase. This gas 
has a short cooling time, and a cooling flow (CF) is formed in the still-forming galaxy. 

In the present paper I assumed that the Mbh — ^fbuige correlation is determined by an 
AGN feedback mechanism that operates during a cooling flow phase at galaxy formation 
(section 4.3). It is determined by the sense that the AGN activity limits the ISM mass that 
is eventually converted to stars by expelling it out of the galaxy. I showed that the Bondi 
accretion cannot operate during that phase (section [2]), and discussed the requirement that 
the ISM be in the hot phase (section EJ). As the radiative cooling time of the hot phase is 
relatively short (eq. [TO)) , to maintain a hot phase the infalling gas should be shocked at a 
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radius of R > 5 kpc. Heating by the AGN facilitates the formation of the hot phase. As 
with CFs in clusters of galaxy, the heating and ejection of the ISM is done by jets, rather 
than by radiation. 

The short radiative cooling time implies the formation of a CF in the the inner region of 
the newly formed galaxy, but one that is substantially heated by the AGN activity. As Bondi 
accretion fails, the feeding of the SMBH is done via two channels. Like in the cold feedback 
mechanism in clusters of galaxies (Pizzolato & Soker 2005, 2010; Soker 2006, 2008b), cold 
blobs are falling from an extended region. The second channel is a cold supersonic inflow 
in the inner < 1 kpc of the galaxy, where cooling time is shorter than the inflow time; such 
an inflow was studied by Soker & Sarazin (1988). The AGN jets can be efficient enough 
to expel a large fraction of the inflowing gas and the gas in the hot phase (Soker 2009) out 
of the galaxy. Some of the cooling gas will form stars and feed the SMBH. A cooling flow 
where a large fraction of the cooling gas is expelled form the inner region (and some fraction 
forms stars, and a small fraction is accreted by the SMBH) is termed a moderate CF model 
(Soker et al. 2001, Soker & David 2003). In section I4T21 the moderate CF model in clusters 
of galaxies and at galaxy formation are compared, and summarized in Table 1. 

The moderate CF model proposed here at galaxy formation can be applicable to the 
sample of obscured AGN studied by Brusa et al. (2009). The median value of -^agn / -^Edd 
in their sample is ~ 2 — 10%. This gives a SBMH growth time of ~ 10 9 yr, similar to the 
growth time of the stellar population due to star formation. I propose that a moderate CF 
exists in these galaxies during the high star formation rate. 

I thank an anonymous referee for helpful comments. This research was supported by 
the Asher Fund for Space Research at the Technion, and the Israel Science foundation. 
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